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Abstract 



A double exchange model for degenerate Cg orbitals with intra- and inter-orbital 
interactions has been studied for the electron doped manganites Ai_a,Ba,Mn03 (x > 
0.5). We show that such a model reproduces the observed phase diagram and orbital 
ordering in the intermediate bandwidth regime and the Jahn-Teller effect, considered 
to be crucial for the region x < 0.5, does not play a major role in this region. Brink 
and Khomskii have already pointed this out and stressed the relevance of the anistropic 
hopping across the degenerate eg orbitals in the infinite Hund's coupling limit. From 
a more realistic calculation with finite Hund's coupling, we show that inclusion of 
interactions stabilizes the C-phase, the antiferromagnetic metallic A-phase moves closer 
to X = 0.5 while the ferromagnetic phase shrinks. This is in agreement with the recent 
observations of Kajimoto et. al. and Akimoto et. al. 

PACS Nos. 75.30.Et, 75.30.Vn 

The perovskite manganites have been at the centre of attention [0, |^ recently as systems 
like Pri_a;Ca2,Mn03, Lai_a;Ca2,Mn03, exhibit colossal magnetoresistance (CMR). These two 
extensively studied systems have relatively low bandwidths. The observation of CMR even 
in the wider bandwidth materials like Pri_a;Sra;MnO3 0] and Ndi-ajSr^jMnOs for x > 0.5 
have prompted a series of careful experimental work on the magnetic phase diagram of all of 
these systems particularly in the region x > 0.5 (the so called electron doped regime). This 
region shows a rich variety of magnetic and orbital ordering and only recently the systematics 
of the phase diagram with externally controlled bandwidth have begun to emerge^. 

In their study of Pri_a;Sra;Mn03, Kajimoto et. al. have summarized the nature of 
magnetic ordering for a series of manganites having different bandwidths (Fig.l in ref. 0). 
They observe that there is no CE phase in many of these systems and the sequence of 
phases in the electron doped region for moderate to large bandwidth systems follows the 
order ferromagnetic (F) — >■ A-type antiferromagnetic (AFM) — C-type AFM and finally to 
G-type AFM phase. The F phase close to x = 0.5 is very narrow and survives for systems 
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with bandwidths above Pro.sSro.sMnOs (e.g., in (Lao.5Ndo.5)i-xSra; MnOs there is a small 
sliver of F metallic phase ). A-phase exists in a small region while the C-phase covers the 
widest region in the phase diagram. The general trend is that with decreasing bandwidth 
the F-phase reduces, A-phase moves closer to a; = 0.5 while the C-phase grows. It is also 
observed that the gradual building of AFM correlations, starting from x = 0.5, is preempted 
by the orbital ordering in the A and C phases]^, There does not seem to be any convincing 
evidence in favour of phase separation in this region[p. 

In the absence of Jahn- Teller (JT) splitting the two e^ orbitals of Mn ion are degenerate. 
The doped manganite Ri-j^A^MnOs has y = 1 — x number of electrons in the Cg orbitals and 
the filling, therefore, is |. In the foregoing, we restrict ourselves to the region x > 0.5 i.e., 
y < 0.5. At the x = 1 end, the band is empty and the physics is governed entirely by the 
exchange between the t2g electrons in the neighbouring sites. On doping, the kinetic energy 
of electrons in the e^ levels begin to compete with the AF superexchange (SE) between 
neighbouring t2g spins via Hund's coupling and this leads to a rich variety of magnetic and 
orbital structures. A model incorporating this physics has recently been proposed by Brink 
and Khomskii]^] (hereinafter referred to as BK). 

The model BK used for the electron-doped manganites contains three terms 

H = JafY. Si-Sj - JhY. Si-^i - ttd4,a,aCj,f3,a (1) 

<ij> i <ij>a,a,l3 

The first term represents the AF exchange between t2g spins, the second term is the 
Hund's coupling between t2g and e^ spins at each site and the third one provides hopping 
between the two orbitals |jlO|] {a, (3 take values 1 and 2 for d^-i-y-^ and d322_^.2 orbitals, corre- 



sponding to the choice = in Ref.[^). 

The existence of the degenerate e^ orbitals with the asymmetric hopping integrals t^^ 
makes (1) very different from the usual DE model|I[]. BK treated the i2g spins quasi- 
classically and Jh was set to infinity. Canting in the x-z plane was included through the 
effective hoppings t^y = tcos{9xy/2) and = tcosiOz/l). Here O^y {6^) is the angle between 
nearest neighbour i2g spins in the x-y plane (z-direction). The superexchange energy per 
state, then, is Ese = ^ C^cosOxy + cosOz). At this level of approximation, the prob- 
lem reduces to solving the 2x2 matrix equation \\ta(3 — e5a/3|| = for a system of spinless 
fermions and minimizing the total energy with respect to 9xy and 9z- In the uncanted state, 
9xy = ^2 = implies F phase, Oxy = 0^ = 71 G-type, 9xy = vr and 9z = C-type and 6xy = 
and 6z = 7T A-type AFM phases. 

Remarkably, the phase diagram obtained by BK based on such simplifying assump- 
tions indeed shows the different magnetic phases seen experimentally in Ndi-^SrajMnOs 
and Pri_a;Sra;Mn03 in the region x > 0.5 although the G-phase, expected for the nearly 
empty band (close to x = 1) and the F state at Jaf were, however, not recovered. 
The limit of infinite Hund's coupling which BK worked with is somewhat unphysical for the 
manganites considered |T2|, P^ , |T^. In a more realistic treatment Pai|T^ considered the limit 



of finite and succeeded in recovering the G and F phases. 

From these results BK argue that the degeneracy of Cg orbitals and the anisotropy of 
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hopping are crucial and the JT effect not quite as relevant since the number of JT centres 
is low in the range of doping considered. The effect of disorder, completely ignored in this 
model, does not seem to play a major role in the magnetic phase diagram[l^. 

Neither of the treatments of BK and Pai include the interactions present in the system, 
namely the inter- and intra-orbital Coulomb interactions as well as the intersite Coulomb 
interaction p!2|, O, O]. Although for low doping the interactions are expected to be ineffec- 



tive, with increase in doping they tend to localize the carriers and preferentially enhance the 
orbital ordering. This affects the F-phase and alters the relative stability of A and C phases. 
It is, therefore, necessary to include them in the Hamiltonian (1) and look for their effects 



on the phase diagram. A natural extension to the model (1) is then ||12|| 



int 



i(T(T' 



(2) 



Here U and U' represent the intra- and inter-orbital Coulomb interaction strengths. For 



the systems concerned, we are not looking for the charge ordered statesflg] and neglect longer 
range interactions. We take the interactions U and U' as parameters as in[|ll|, |n[ while in 
reality, these are related through the Racah parameters (see [0 and references therein). We 
treat the spin system quasi-classically, but unlike BK we work at finite Hund's coupling. 
In the uncanted states, we assume for the t2g spin Si = So exp(zQ.ri), where the choice of 
Q determines different spin arrangements for the core spins. In the infinite J// limit, the 
Qg electron spins would be forced to follow the core spins leading to the freezing of their 
spin degrees of freedom. A finite value for ^h-, however, allows for fluctuations and the spin 
degrees of freedom, along with anisotropic hopping across the two orbitals, play a crucial 
role. 

Let us first look at the Hamiltonian (1) i.e., set U = U' = 0. For t2g spin configurations 
described above, it reduces to 



H 



'^kQ(T'^k/3(7 



k,a,/3,a- 



JhSq C^a|Ck+Qat + JhSq ^ cj^Q,|Ck+Q. 
k 



(3) 



k,a 



where we have followed the notation in|TD[ for e'^^ . 

We calculate the ground state energy by diagonalization of the above Hamiltonian in a 
finite momentum grid (numerical results converged by a grid size 64x64 x 64) as a function of 
Jh for a range of values of Jaf- The magnetic phase diagram for J^F'S'o = 0.05 in the Jr — x 
plane is shown in Fig. la. The phase diagram in J^f — x plane for J/^Sq =10 is plotted in 
Fig. 2a. All energies are measured in units of the hopping t. There is no general agreement on 

From photoemission and optical studies [p 



the values of the parameters involved|12 



1311 and 



1.5 - 2 



LDA analysis JTSP one can glean a range of typical values Q.leV < t < 0.3el^, J// ^ 
eV and J^^ — 0.03t — O.lt (Maezono et. al.|jl9[ quote a lesser value of J^^ = O.Olt). 

At the X = 1 end, with empty e^ orbitals, the only contribution to energy comes from 
the SE interaction leading to the G-type AFM phase. On doping by electrons the C-phase 
appears first with orbital ordering (of the orbitals) along the z-direction. The stability of 
the A-phase comes from the ordering of dx^-y^ orbitals in the xy plane. The gain in kinetic 
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energy due to the planar and one dimensional orbital order, induced by the anisotropic 
hopping integral, more than offsets the loss of SE energy. The orbital order drives the 
corresponding magnetic order as well - in the A-phase the spins have planar FM order and 
AFM order along the z-direction whereas in the C-phase it is reversed. The 3D magnetically 
ordered F- and G-phases show no orbital ordering. The phase transitions are therefore 
characterised by the density of states (DOS) reflecting the underlying 1, 2 and 3 dimensional 
characters of the different phases. The phase diagrams are shown in Figs, la and 2a, for 
typical values of the parameters JhSq = 10 and JafSq = 0.05. The sequence of G, C, A and 
finally the F phase with complete alignment of spins is observed. 

In the infinite Hund's coupling limit the "wrong" spin sector of the Hilbert space was 
projected out by BK. When the system is doped, spin canting is the only channel for de- 
localization of doped carriers, albeit with a loss in SE energy. At finite Jh, however, the 
wrong spins are no longer as "costly" and canting is expected to reduce. Canting is included 
through the choice Si = So(sm6'j, 0, cos^j). The Hund's coupling term in the Hamiltonian 

becomes Hh^nd = -JHSoT,i,aCos6i{cl^^Ciai - c|„|Cia|) - J//5'o s^^^i(cLT'^i"i + ^al^ia^)- 
In this case the different magnetic phases need to be defined at the outset. The convention 
used by BK to define the magnetic phases are: it is A-type when O^y < 9^ and C-type when 
dxy > 9z- In the canted G and F phases O^y and 9^ are close to 180° and 0° respectively, 
although, it is obvious that the canted G-phase and A-phase are synonymous in a certain 



region. However, orbital order can be used to delineate the two phases [pi|. 

Proceeding as before, the ground state energy for different 9xy and 9^ is obtained. The 
qualitative phase diagram is very similar to the uncanted case except for small shifts in the 
phase boundaries (the shifts are small unless Jh is large) and agrees with Pai|T^. We 
show in Fig.3 the angle of canting for both 9z and 9xy deep inside the G-phase at x = 0.98. 
The angles in Fig.3 represent deviation from 180°. There is hardly any canting in 9^ while 
in 9xy, there is no significant canting for low and it is about 10° only for large J^- We 
note that experiments P, have so far not been able to detect any significant canting in A- 
and C-phases. Even in the G-phase, certain systems appear to show little canting. 

The interactions (both intra- and inter-orbital) are treated in the mean-field (MF) theory 
and a self-consistent calculation has been performed. Self-consistency is achieved when all 
the averages < rii,o-,a > and the ground state energy converge to within 0.01%. Figs. la, b 
and 2a,b show the modifications in the magnetic phase diagram by the inter-orbital {U') 
interaction in the x — Jh plane at JafSq = 0.05 and in the x — Jaf plane for JhSq = 10. 
Although we obtained the phase diagram for several values of U', in Figs. 1,2 we only give 
representative ones for demonstration. On increasing U', the F-phase starts shrinking fast, 
the C-phase gains in size while the G-phase remains almost unaltered. This is primarily 
because of the enhanced orbital ordering in the A- and C-phases driven by the inter-orbital 
repulsion and the low dimensional nature of the DOS in these phases. As discussed earlier 
the AFM A and C phases are driven by orbital ordering and in the presence of U', the one 
dimensional order leading to the AF instabilty in the C-phase grows faster. Close to the 
X = 1 end the electron density is very low, there are almost no sites with both the orbitals 
occupied and U' is therefore ineffective. At the other end, however, the density is higher and 



4 



the F phase has preferential occupation of one spin at both the orbitals. Hence this phase 
is affected drastically by the inter-orbital repulsion. 

It is known [11, 12 1 that at the level of MF theory the intra- orbital repulsion U between 
opposite spins mimics the effect of Jh- As we are working with quite low densities (actual 
filling < 0.125), and the relevant Jh values being large, we find almost no observable effect 
of U on the phase diagram (except for very low Jh where again the changes are small). 

In the absence of interactions there is orbital ordering in both A- and C-phases. The 
presence of U' enhances this ordering. We calculate the orbital densities in both A- and 
C-phases in the respective ground states and show the results in Fig. 4. In the figure, we 
have plotted actual orbital occupancies (the sum of the occupancies of the two orbitals will 
be ^^) for different U'. It is evident from the figure that in the A-phase the dx2-y2 orbitals 
are predominantly occupied, while in the C-phase the d22 orbitals have higher occupancy. 
As U' increases, the orbital ordering is enhanced. This is shown in Fig. 4 for three values of 
U' in the A and C phases in their respective regions of stability as a function of doping. Note 
that as X increases (density of electron decreases), the effect of U' on the orbital occupancies 
becomes less pronounced and the curves for different U' merge as expected. We also show 
the orbital occupancies as a function of U' in Fig. 5 in the regions where A- and C-phases are 
stable and the effect of U' is noticeable in both the A- and C-phases. The orbital densities 
in C-phase attain their saturation values by [/' ~ 8. Since we are interested in the region 
X > 0.5, we have not plotted orbital densities in A-phase beyond U' = 8 - above this value 
A-phase shifts below x = 0.5 at J// = 5 (see Fig.l). 

The present calculations produce results that agree with BK and Pai for U=U'=0. We 
are able to recover the G-phase at x ~ 1 and we also obtained the F phase for Jaf — 0. On 
inclusion of inter- and intra-orbital interactions, the topology of the phase diagram remains 
the same. The C-phase grows at the expense of F phase while the G-phase remains unaffected 
with increasing This scenario is borne out in the bandwidth controlled experiments of 
Akimoto et. al.p| and the schematic phase diagram obtained by Kajimoto et. al|P]. Our 
results qualitatively agree with the earlier work of Maezono et. al.[|19| as well. They included 
correlations in an MF treatement, but did not get the A-phase close to a; ~ 0.5 observed 
experimentally. Although Fig. 15 in Maezono et. al. |^ resembles (with a vanishing A- 
phase close to a; = 0.5) our Fig. 2, a comparison is quite difficult owing to the very different 
choice of the parameters (it is also not possible to separate the effects of Coulomb and 



exchange interactions in their work). In their Monte Carlo treatment, Shen and Tingp2 



considered an effective model and obtained a phase diagram. However, the C-phase in the 
region 0.6 < x < 0.9 does not come out of their work. Hotta et. al.[ll] have compared 
exact diagonalization results in one dimension with MF theory and found the agreement 
to be good. Our MF calculations also suggest that the qualitative trends obtained are in 
good agreement with the physically expected and experimentally observed results in the 
manganites. 

We note that the value of U' for which the F-phase disappears from the region x > 0.5 in 
our calculation is about ^ = 12. Depending on the value of t, corresponding U' is between 
1.8 - 3.6 eV. This value is somewhat on the lower side for the range of values available in 
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literature (the range varies between 3-10 eV)|T^, 0, 0]. Although the available values are 
actually bare values and in phases like F and A, they are bound to go down owing to metallic 
screening (a treatment of which is beyond the scope of this work) - a problem faced in all 
theories of correlated systems across a metal-insulator transition. In the foregoing, we have 
assumed that increase in interactions is qualitatively equivalent to reduction of bandwidth, 
while in reality, the interactions play more complex roles in addition to charge localization 
which are not included in our calculation. 

In conclusion, we have included orbital correlations in a degenerate double exchange 
model proposed by Brink and Khomskii for the electron doped, intermediate bandwidth 
manganites. We observe from a generalized mean-field calculation that the phase diagram 
captures most of the qualitative features seen experimentally. The orbital orderings ob- 
tained are in good agreement with experimental observations. It also agrees with the trends 
observed across several manganites with changing bandwidths. It would be interesting to 
include JT coupling and extended range Coulomb terms in the model and observe their 
effects particularly close to the x = 0.5 region. 
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Figure captions 



Fig. 1. Magnetic phase diagram in doping (x) - Jh plane with (a) U' — and (b) U' — 8.0. 
All energies are measured in units of t. 

Fig. 2. Magnetic phase diagram in doping (x) - Jaf plane with (a) U' = and (b) U' = 8.0. 

Fig. 3. Canting of the angles 9xy and 9z in degrees as a function of Jh {JafSq — 0.05). 

Fig. 4. Orbital densities as a function of doping x for three values of U' — 0, 4, 8. The filled 
symbols are for (i^2 and open symbols for d^2_y2 orbitals. The vertical dotted lines 
represent the boundary between A- and C-phases for different U'. We choose JhSo = 5 
here in order to have stable A- and C-phases for a reasonable range of x (see Fig.l) for 
all three U' values. JafSq was kept at 0.05. 

Fig. 5. Orbital density versus U' in (a) A-phase at a; = 0.5 and (b) C-phase at x = 0.65. The 
dotted lines are for d^^ and solid lines are for dx^-y^ orbitals. JhSq and JafS^ were 
same as in Fig. 4. 
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